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We study three tight-binding models of two-dimensional fermionic optical lattices. We find that 
these systems have a novel composite symmetry, i.e. gauge-color-translation (GCT) symmetry, 
which consists of a fixed translation transformation, a color transformation and a fixed local U(l) 
gauge transformation. The corresponding symmetry operator is antiunitary. We predict that, due to 
the protection of GCT symmetry, the conduction and valence bands are degenerate at some isolated 
momenta. These degenerate points form Weyl nodes in the Brillouin zone. The quasiparticles and 
quasiholes near the Weyl nodes can be considered as Weyl fermions and have chirality characterized 
by a winding number. With the protection of GCT symmetry and the winding number, the Weyl 
nodes occur in a very large range of parameters. 
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Introduction. — The research on topological phases is 
becoming an increasingly important theme in condensed 
matters. Topological matters are classified according to 
topological invariants rather than symmetries [HQ. More 
recently, physicists have found that, besides the gapped 
systems, gapless systems can also support topological 
phases, i.e., topological semimetals [3l4l 2] . Semimetals 
are characterized by the presence of point nodes in the 
band structures, that is, two distinct bands touch each 
other at some isolated points in momentum space. 

Weyl semimetal is a class of topological semimetals. 
For Weyl semimetals, near the band-degeneracy points 
(also called Weyl nodes), the single particle Hamiltonian 
can be written as h(k) = Y^ij VijkiCTj, with k and a be- 
ing the wave vector and the Pauli matrix, respectively. 
That is to say, near these touching points the dispersion 
relation is linear and can be described by a massless two- 
component Weyl Hamiltonian. Near the Weyl nodes, the 
quasiparticles are Weyl fermions and have chirality char- 
acterized by a winding number denned as 



= sgn[det(^)] = ±1. 



(1) 



Due to the theorem of Nielsen and Ninomiya, Weyl nodes 
always occur in pairs of opposite chirality 13]. In Weyl 
semimetals, Weyl nodes with opposite chirality are sepa- 
rated from each other in momentum space and thus can 
not be destroyed by mutual annihilation of pairs with 
opposite topological charges. 

In three dimensions, the band degeneracy at Weyl 
nodes can be accidental Q or protected by symmetries. 
The robustness of the accidental degeneracy depends on 
its co- dimension (EJ. To achieve a two-band accidental de- 
generacy, three real parameters are required to be tuned. 
Thus, such accidental band degeneracies are vanishingly 
improbable in one and two dimensions if there are not 
additional symmetry constraints [15j. Therefore, in two 
dimensions, the band degeneracy at Weyl nodes must be 
protected by symmetries. In general, the band degen- 



eracies are protected by point groups or time reversal 
symmetry. In this Letter, we will show that, besides 
point groups and time-reversal symmetry, there is an- 
other kind of symmetry, gauge-color-translation (GCT) 
symmetry, to protect band degeneracies in condensed 
matters. GCT symmetry protected Weyl nodes are ubiq- 
uitous in condensed matters, but which are rarely stud- 
ied by physicists before. GCT symmetry is a composite 
symmetry; its symmetry operator is a combined opera- 
tion and consists of a translation transformation, a color 
transformation and a local gauge transformation in se- 
quence. Sometimes, the local gauge transformation is 
not necessary, which can be regarded as a trivial case 
with a global gauge transformation. GCT operator has 
the form as T = GCT T1 where T T is a translation opera- 
tor such as T T r = r + r, C is a color transformation like 
cr x K with K and cr x being the complex conjugation and 
pauli matrix, respectively, and G is a fixed local gauge 
transformation. Here, we stress that GCT symmetry is 
an antiunitary symmetry, i.e., its symmetry operator is 
antiunitary, because it includes the complex conjugation 
K. GCT symmetry operator Y and the space group of 
the lattice can make up a generalized color space group, 
where a color transformation is always followed by a lo- 
cal gauge transformation. However, only GCT symmetry 
is responsible for protecting Weyl nodes. Thus, here we 
mainly consider GCT symmetry. 

In order to present GCT symmetry protected Weyl 
nodes, we first consider a two-dimensional rectangular 
optical lattice as shown in FigOJa), denoted as Model 1. 
We will show that Weyl semimetals appear in this model. 
In order to show ubiquity of GCT symmetry protected 
Weyl nodes in condensed matters, we will also briefly 
investigate and discuss two other two-dimensional rect- 
angular lattices, denoted as Models 2 and 3, respectively, 
as shown in FigQJb) and (c), which also support the ex- 
istence of Weyl semimetals. These models can experi- 
mentally be realized and detected in optical lattices with 
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FIG. 1: (Color online). Schematic of the rectangular optical 
lattices and the designed phase factor (denoted by arrows) for 
(a) Model 1, (b) Model 2, and (c) Model 3. Here, the green 
and blue balls represent the lattice sites of sublattices A and 
B, respectively. 



laser-assisted tunneling as proposed in Reference [16j. 

Model. — Model 1 consists of two sublattices denoted 
as A and £?, respectively. The two sublattices have the 
lattice spacings 2d x and d y in the x and y directions, 
respectively. For sublattice along y direction, there 
exists a phase factor between two neighbor lattice sites. 
For each sublattice, the primitive lattice vectors are de- 
fined as a\ = (2^,0) and a 2 — (0,d y ). The lattice 
vectors can be expressed as R n = n\a\ + n 2 a 2 with rt\ 
and n 2 being integers, we chose a lattice site in sublat- 
tice A as the origin, then the lattice vectors in sublattice 
A and sublattice B can be rewritten as = R n and 
R^ = R n + ai/2. 

The tight-binding Hamiltonian for Model 1 can be 
written as, 

+ Ji,i+xa\bi+x + Ji+x,i+2xb\j r3 ,Cli+2x + H.c], (2) 



+y = ty for i G A and t y e 11 for 



with Ji,i+ X — t x , J 2 
i G B. Here, is the annihilation operator that destructs 
a particle with the Wannier state wf(r) = w^(r — Rf) 
located at the site i in sublattice A and bj the annihilation 
operator destructing a particle with the Wannier state 
w f( r ) = w o( r — Rf) located at the site j in sublattice 
B; the subscript i = (i x ,i y ) is the coordinate for the sites 
of sublattice A with the primitive lattice vectors a\ and 
a 2 as the bases; x = (1/2, 0) and y = (0, 1) represent the 
vectors between the neighbor lattice sites in the whole 
lattice; < 7 < 2ir is a phase factor along with tunneling; 
t x and t y are the hopping amplitudes along the x and y 
directions, respectively. 

This system is invariant under a GCT transformation, 
i.e. H = Ti^r -1 , where the symmetry operator T is 



T = (e^yya x KT ai/2 , 



(3) 



where (e ll ) ly is a local U(l) gauge transforma- 
tion. The corresponding inverse operator is T _1 = 
(e^v a x KT~* /2 . The GCT operator has the character 
as follows, 



1 a 1 5 



(4) 



where T ai = {E\ ai }. 

Band structure and Block functions. — In the recipro- 
cal lattice, the reciprocal lattice vectors are defined as 
K rn = mi 61 + m2 62, where mi and 777,2 are integers, and 
&i = (tt/4, 0) and 62 = (0, 2k /d y ) are the corresponding 
primitive reciprocal lattice vectors. Thus, the Brillouin 
zone is —7r/2d x < k x < ir/2d x , —7r/d y < k y < 7r/d y . The 
corresponding Bloch functions for the A and B sublat- 
tices can be written as ip^ (r) 

and Tjj*(r) = ^- e 
respectively, where N is the number of lattice sites in 
each sublattice. Correspondingly, we have the annihi- 
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lation operators as a& — 

~ lk ' R i 5 which destruct a particle in the Bloch 
states an d V'fc respectively. We define the 

two-component annihilation operator as 77^ = :^[ a fc + 
ifrfc, za/g + bk] T • The Hamiltonian (|2j) can be rewritten as 
H = Z k vln(k)r lk with 

'H(fc) = — t 2y [cos(/c 2y (i 2y ) + cos(kyd y — 7)]/ — 2^ cos(k x d x )a x 

— ty[COS(kydy) ~ COS(kydy ~ J ) ] Gy , (5) 

and cr x and are the Pauli matrices and / the unit 
matrix. 

Diagonalizing Eq.(j5]), we obtain the dispersion re- 
lation as E(k) = —ty[cos(kyd y ) + cos(k y d y — 7)] ± 
t y ^4/ 2 cos 2 (k x d x ) + [cos(/c 2y (i 2y ) — cos(k y d y — j)} 2 with 
/ = tp/t^. Fig E] shows this dispersion relation for the 
cases with 7 = tt and 7 = 3tt/5. We assume that 
[c^,c^] T is the eigenvector of the Hamiltonian (|5j) cor- 
responding to eigenenergy E(k), then the Bloch function 
of the whole lattice can be expressed as 



ui(r) 
u 2 (r) 



Ahr 



CfeVfe (r) - 4 2) i^(r) 
-4 1) #f(r-) + 4 2 Vf(r-) 



where ^(r) = i^(r + i£ n ). In the momentum space, the 
Bloch state ^fe(r) and eigenenergy 22(fc) are periodic for 
reciprocal lattice vectors, i.e. ^fe(r) = ^ r fc + K: rn (r) and 
E(k) = E(k + K m ). 

Weyl semimetals. — From the dispersion relation, we 
find that the conduction and valence bands are touched 
at Wi = (-7T/2d x ^/2d y ), W 2 = (7r/2d x ,~f/2d y ), W 3 = 
(—7v/2d x ,j/2dy — 7v/d y ) and W4 = (ir/2d x , j/2d y — 7r/d y ) 
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FIG. 2: (Color online). The energy dispersions of Model 1 for 
(a) 7 = 7r and (b) 7 = 3ty/5. Here, = t y is assumed. 
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in the Brillouin zone. Near these degenerate points, the 
single-particle Hamiltonian ([5]) can be linearized as 

h(p) = ±2t y cos(j/2)I±2tyfd x p x a x 

±2ty Sm(~i/2)dyPy(Jy, (7) 

where the signs before the three terms make up four com- 
binations such as (—,—,+), (—,+,+), (+,—,—), and 
(+,+,—) representing the linearized Hamiltonian near 
the touched points Wi, W2, W3, and W4, respectively. 
The Hamiltonian ([7]) has the form of Weyl's Hamilto- 
nian and the four Weyl nodes which are located at the 
boundary of the Brillouin zone. Because the relations 
W2 = W\ + 61 and W4 = Ws + 61, W\ is equiva- 
lent with W2, and W3 equivalent with W4, i.e., there 
are only two unequivalent Weyl nodes in the Brillouin 
zone. From the the definition of the chirality as Eq.(pQ), 
we have w = — 1, +1, — 1, +1 for Weyl nodes Wi, W2, 
W3, and W4, respectively. We find that W\ and W2, 
W3 and W4 are equivalent points but have opposite chi- 
rality. Thus, we conclude that the chirality is not abso- 
lute in our model, which can be changed through a gauge 
transformation. However, due to the theorem of Nielsen 
and Ninomiya[l3|, in any gauge, the two unequivalent 
Weyl nodes have opposite chirality. 

The Weyl nodes can also interpreted as topological de- 
fects of the planar vector field in momentum space like 
h = (-2t y f cos(k x d x ), -t y [cos(kyd y ) - cos(kyd y - 7)]), 
which is denned by using the coefficients of the Pauli ma- 
trices in Eq. (j5j) . The corresponding topological invariant 
is the winding number denned as w = § c |^ • F(k) with 

F(k) = ft V (ft) - ft V (ft) 0. These vortices have 
the values of winding number ±1, which is consistent 
with the chirality denned in Eq.(pQ). Fig |3] shows the vec- 
tor fields F(k) for the cases with 7 = tt and 7 = 37r/5, 
which manifest the vortex structures indicating the ap- 
pearing of Weyl nodes. 

Symmetry protection. — The Weyl nodes in our systems 
manifest strong robustness. Firstly, the Weyl nodes ap- 
pear and the band structure keep gapless for all ranges 
of < 7 < 2tt. Secondly, any values of the lattice con- 
stants d x and d y can support Weyl nodes. Besides the 
protection of topological invariant, these Weyl nodes are 
protected by GCT symmetry, the operator T of which is 



FIG. 3: (Color online). The vortex structures of the vector 
fields F(k) for Model 1 with (a) 7 = tt and (b) 7 = 3tt/5. 
Here, the region surrounded by red lines is the first Bril- 
louin zone; the green balls mark the GCT-invariant points 
in the first Brillouin zone; the vortices indicate the appearing 
of Weyl nodes. 



expressed in Eq.©. 

The GCT operator T acts on the Bloch function (|6|) as 
follows 



u^ > (r-a 1 /2)e ikxdx 
u£ ) *(r-a 1 /2)e ikxdx 



e -i[k x x+(k y -j/dy)y] 

(8) 



where x = i x d x and y = i y d y . Because T is the sym- 
metry operator of the system, ty' k ,(r) must be a Bloch 
wave function of the system. Thus we obtain u^,\r) = 
u^ > {r-a l /2)e ik ^ and v${r) = (r-a 1 /2)e ikxd ^ 
kg, —— k x and ky 7 / dy 

From Eqs. (j6]) and (j8j), it is easy to show that the 
operator T has the effect when acting on the wave vec- 
tor k = {k x ,ky) as T : {k XJ ky) > (—k x ^jdy — ky). 
If k' = k + K mi then we can say that k is an in- 
variant point under GCT transformation. In the Bril- 
louin zone, the GCT-invariant wave vectors are — 
(±7r/2d xn /2dy),G 3A = (±7r/2d Xl ( 1 /2-7r)/dy), G 5 = 
(0, *y/2d y ), Gq = (0, (7/2 — 7r)/d y ), which are marked by 
green balls in Fig 131 For the case k' = k + K m , i.e. at the 
GCT-invariant point G;, we have T^ Gt (r) = Sfr' Gi (r). 
Thus, (?) and ^' G .(r) are both the eigenstates of 
Hamiltonian H and have the same eigenenergy E(Gi). 
Considering Eq.©, we have T 2 ^/ Gi (r) = T ai ^ Gt (r) = 
e~ 2lGixdx ty G .(r), where G{ x is the x component of G^. 
We define (ip, ip) as the inner product of the two wave 
functions ip and cp. The antiunitary operator T has the 
property that (Tip^Tcp) = (-0, = (<p,ip). Therefore, 
we have 

(*' Gl ,* G <) = (r^,r^ G j = (r^,r 2 ^) 

= e- 2iG ^ dx (^' G .^ G .). (9) 

For the case with Gi X = 0, such as G5 and Gq : 
(^ G . : ^gJ is unconstrained. For the case with d x ^ 0, 
such as ^1,^2,(^3 and G4, we obtain the solution 
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(^' G . , ^Gi) = 0, that is to say, ^' G . and are orthogo- 
nal to each other. Therefore, we arrive at the conclusion 
that the system must be degenerate at points = 
(±7r/2d x ,j/2d y ) and G 3A = (±ir/2d x , (7/2 - n)/d y ) in 
the Brillouin zone, which are just the Weyl nodes Wi, 
W 2 , W 3 , and W 4 . 

Other models — Besides the above model, there are 
other models, such as Models 2 and 3 shown in FigHJ that 
have similar composite symmetry and the corresponding 
symmetry protected Weyl nodes. Model 2 has the same 
sublattices and the same Brillouin zone with Model 1. 
Thus, its tight-binding Hamiltonian is similar with Eq.(|2|) 
but with Ji,i+ X = t x e~ 1 ^ , Ji,i+ y = t y e~ %1 for i G A and 
t y e 11 for i G 5, where < /?, 7 < tt. This model has the 
spectrum as E'(k) = —t y [cos(k y d y —j) + cos(k y d y -\-j)]± 
tyV^f 2 ^os 2 (k x d x - 0) + [cosp^ - 7) - cos(k y d y + 7)] 2 
with / = t x /ty. The Brillouin zone is —n/2d x < k x < 
7r/2d X: —7Y/d y < k y < n/d y . Weyl nodes occur at 
(-7r/2d x + /3/d x ,0) and (-n/2d x + (i/d x ,±7r), among 
which there are only two unequivalent points. This model 
has a composite symmetry as r' = (e 2l ^Y x a x KT ai /2 and 
p/2 _ e 2*/5j> a ^ j n momentum space, the wave vectors are 
transformed as follows T' : (fc^, fc y ) — (2/3 jd x — k x , —k y ). 
For the < f3 < tt/2 case, the r'-invariant points in the 
Brillouin zone are (— tt/2 + /3/d x , 0), (—tt/2 + f3/d x , ±7r), 
(f3/d x , 0), ((3/d x , ±7r) and for the tt/2 < f3 < n case they 
are(-7r/2 + /3/^,0), (-n/2 + P/d x ,±ir), ((P-n)/d x , 0), 
((/? — n)/d x , ±7r). Here, only at points (—tt/2 + P/d x , 0), 
(— 7r/2 + (3/d x ,±7r), the Bloch wave functions \I/fc(r) 
and T / ^fc(r) are orthogonal, i.e, the degeneracies at 
these points are protected by GCT symmetry, which 
is consistent with the observation of the dispersion of 
Model 2. 

Model 3 also have two sublattices and each sub- 
lattice has the unit lattice vectors a'/ = (d Xl —d y ) 
and a 2 ; = (d x ,d y ). In reciprocal lattice space, the 
unit wave vectors are b'{ = (ir / d x , —ir / d y ) and b f 2 ' = 
(ir / d x , ir / d y ) . The tight-binding Hamiltonian can be 
written as H" = -J2ieAl J h^x a l b i+& + Ji,i+ya\b i+y + 

+ ^i+£,i+2x^ + ^^i+2x + H.C.] with 

Ji,i+ X = t x e~ %1 for z G A and t x e %1 for z G 5, 
Ji 5 i+y = ^e* 7 for i E A and t y e~ %1 for z G 5, where 
< 7 < 7r/2. The spectrum has the form as E" (k) = 

±2<^cos 2 7^5j_ + sin 2 7^ 2 _ with Q + = t x cos(k x d x ) + 
^ cos(kyd y ) and ^_ = ^ cos(k x d x ) — t y cos(k y d y ). Weyl 
nodes occur at (±7r/2d x ,±7r/2d y ) in the Brillouin zone, 
among which there are only two unequivalent Weyl 
nodes. The symmetry operator is V" = o- x KT( a »+ a »y 2 
and its square is V" 2 = T a // +a //. In momentum space, the 
wave vectors are transformed as follows, r" : (k x , k y ) — >• 
(—k x ,—ky). The T^-invariant points in the Brillouin zone 
are (±ir/2d x ,±ir/2d y ) and (0,0). But only degeneracies 
at (±7r/2d x ,±7v/2d y ) are protected by GCT symmetry, 
which is consistent with the dispersion calculated. 



Remarks and conclusion. — The models investigated 
above can be realized in optical lattice with laser-assisted 
tunneling [l6|. The GCT symmetry is related with vari- 
ous branches of physics. For example, besides in optical 
lattices, Model 2 with /3 = 7 = 7r/2, = d y and Model 
3 with j = 7r/4,d x = d y are just the 7r-flux states in spin 
liquid 1 17j and high temperature superconductors [18| , re- 
spectively. Three-dimensional Weyl semimetals investi- 
gated in Reference [10[ are also protected by GCT sym- 
metry. Thus, GCT symmety are ubiquitous in physics, 
in particular condensed matter physics, and has a close 
relation with the existence of Weyl nodes. 

In summary, we have studied the three optical lattice 
models that support the existence of Weyl semimetals. 
We found that there exists a novel symmetry, GCT sym- 
metry in these models. GCT symmetry is an antiuni- 
tary symmetry and responsible for point degeneracies of 
Weyl nodes. We have discussed the dispersion relations 
of these models and analyzed the chirality of Weyl nodes. 
Our study on the relation between GCT symmetry and 
Weyl nodes gives a clue to search for Weyl semimetals 
and an insight on interpreting Weyl semimetals. The 
study of GCT symmetry is also helpful to understanding 
the essence of spin liquids and high temperature super- 
conductors. 
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